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Abstract 

We study the late time evolution of positively curved FRW models 
with a scalar field which arises in the conformal frame of the R + 
aR 2 theory. The resulted three-dimensional dynamical system has 
two equilibrium solutions corresponding to a de Sitter space and an 
ever expanding closed universe. We analyze the structure of the first 
equilibrium with the methods of the center manifold theory and, for 
the second equilibrium we apply the normal form theory to obtain a 
simplified system, which we analyze with special phase plane methods. 
It is shown that an initially expanding closed FRW spacetime avoids 
recollapse. 

1 Introduction 



A central question in relativistic cosmology is that of deciding about the past 
and future asymptotic states of cosmological models [T]. However, general 
relativity leads to singularities in the spacetimes of all known cosmological 
models with ordinary matter. Higher order curvature corrections in the grav- 
itational action may rectify the problem and lead to cosmological models free 
from such pathologies, at the cost of diverging from a FRW behavior at late 
times j2]. There is a resurgence of interest in such theories which naturally 
arise in string-theoretic considerations (cf. brane models with Gauss-Bonett 
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terms [21 IU 03 An interesting feature of higher order theories is that 
inflation emerges in these theories in a most direct way. In one of the first 
inflationary models, proposed in 1980 by Starobinsky [7j, inflation is due to 
the R 2 correction term in a gravitational Lagrangian L = R + aR 2 where 
a is a constant. The dynamics of higher order cosmologies is closely related 
to scalar-field cosmologies in general relativity because of their conformal 
equivalence (HI El- There are certain limitations to this procedure related 
to the issue of physical reality of the two metrics involved [TU] and to the 
fact that the conformal transformation may fail to be regular at all points of 
the spacetime. Nevertheless, it is practically useful and investigations in the 
conformal frame have given some interesting results, e.g. the cosmic no-hair 
theorem in quadratic cosmologies [TT| I12|. 

Most of the studies of scalar-field cosmologies with the dynamical sys- 
tems methods are restricted to FRW models (see for example [T51 HI] and 
references therein), although there are important investigations in spatially 
homogeneous Bianchi cosmologies with an exponential potential . In par- 
ticular, for flat FRW models with a scalar field there are some general results 
which do not rely on the particular properties of the potential |Tol IT7] . How- 
ever, the situation is more delicate for positively curved FRW models with a 
scalar field having a potential with a zero local minimum. The main problem 
which confronts us is the following: Can a closed universe filled with ordinary 
matter and a scalar field avoid recollapse? 

In this paper we investigate the evolution of positively curved FRW mod- 
els with a scalar field having the potential which arises in the conformal frame 
of the R + aR 2 theory in vacuum [HJ [nj. The motivation for this choice is 
presented in Section III. The dimension of the dynamical systems involved in 
such models is greater than two and the usual methods of phase plane analy- 
sis cannot be applied. In particular, for nonhyperbolic equilibrium points the 
linearization theorem does not yield any information about the stability of 
the equilibria and therefore, more powerful methods are needed. The center 
manifold theorem shows that the qualitative behavior in a neighborhood of a 
nonhyperbolic equilibrium point q is determined by its behavior on the cen- 
ter manifold near q . Since the dimension of the center manifold is generally 
smaller than the dimension of the dynamical system, this greatly simpli- 
fies the problem. The other general method for simplifying the dynamical 
system is the normal form theory, which consists in a nonlinear coordinate 
transformation that allows to simplify the nonlinear part of the system. Both 
methods are used in Sections III and IV respectively. 
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The plan of the paper is as follows. In the next Section we write down 
the field equations assuming an arbitrary potential, as a constrained four- 
dimensional dynamical system. In Section III, we use the constraint equation 
to reduce the dimension of the system to three and after a suitable change of 
variables we find the equilibrium points of the system and discuss the phys- 
ical meaning of these particular solutions. In particular we show using the 
methods of the center manifold theory that, the equilibrium corresponding 
to the de Sitter solution is asymptotically unstable. In Section IV we find 
the so-called normal form of the dynamical system, which greatly simplifies 
the problem, since two of the equations decouple. We study the qualitative 
behavior of the resulted two-dimensional system and analyze the late time 
evolution of the model. We show that an initially expanding universe avoids 
recollapse. In Section V we apply the same techniques to flat FRW spaces 
filled with a barotropic fluid in the conformal frame of the R + aR 2 theory 
and study the detailed evolution of the models. 

2 Scalar-field cosmologies 

In general relativity the evolution of FRW models with a scalar field (ordinary 
matter is described by a perfect fluid with energy density p and pressure p) 
are governed by the Friedmann equation, 




(1) 



the Raychaudhuri equation, 





the equation of motion of the scalar field, 



<!> + 3-<j> + V' (0) = 



(3) 



and the conservation equation, 



p + 3(p + p) - = 0. 



(4) 



a 
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We adopt the metric and curvature conventions of Ref. [Tj. Here, a (t) is 
the scale factor, an overdot denotes differentiation with respect to time t and 
units have been chosen so that c = 1 = 8nG. 

From Eqs. (HJ-® we see that the state (a, a, p, 0, G IR 5 of the system 

lies on the hypersurface defined by the constraint, (fl|). and the remaining 
field equations can be written as a five- dimensional dynamical system. In 
vacuum, p = 0, the dimension of the dynamical system reduces to four. 

In the literature of scalar-field cosmologies the exponential potential func- 
tion viz. V {(f)) = Voe _A< ^, is the most popular not only because of the variety 
of alternative theories of gravity which predict exponential potentials, but 
also due to the fact that this potential has the nice property that V oc V 
which allows the introduction of normalized variables according to the for- 
malism of Wainwright et al pQ. In flat, k = FRW models for example 
with a scalar field having an exponential potential, introducing the variables 
x ~ 4>/H, y ~ W /H and the time coordinate r = ln(a/ao), enables the 
evolution equations to be written as a two-dimensional dynamical system (cf. 
[18J) and in more general homogeneous cosmologies associated with a scalar 
field, the dimension of the dynamical system reduces by one if the potential 
function is exponential. 

If we set = y, a/a = H the evolution equations (0)-® in vacuum 
become 

a = Ha 

<P = y, 

y = -3Hy-V (0) , 

H = - l -y 2 + k/a\ (5) 

subject to the constraint 

3H 2 + 3k/a 2 = \y 2 + V ((f)) . (6) 
Therefore, the phase space of the dynamical system (JHJ) is the set 

| (a, 0, y, H)eR*: 3H 2 + 3k /a 2 = l -y 2 + V (0) } . 
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3 Curved FRW in R + aR 2 theory: Equilibria 



In the remainder of the paper we assume that the potential function of the 
scalar field is 

yfo) = V„(l-e-V^) a (7) 

which arises in the conformal frame of the R + aR 2 theory [HI EJ- This 
potential has a long and flat plateau. For large values of 0, the potential, V, 
is almost constant, = lim^ +00 V (0) , thus V has the general properties 
for inflation to commence. In J2| it was proved a cosmic no-hair theorem, i.e. 
Bianchi models with ordinary matter satisfying the strong energy condition 
and a scalar field with potential (J7|, asymptotically isotopize. According to 
this picture, the universe started in a homogeneous state and during inflation 
it had enough time to isotropize. 

In order to reduce the dimension of the dynamical system (jSJ) we use the 
constraint (jBJ) to eliminate a. The evolution equations become 

i = y, 

y = -3Hy-V (</>), 

H = -H 2 - l -y 2 + l -V{ct ) ). (8) 

Linearization of (JBJ) near the equilibrium point (0, 0, 0) shows that the Jaco- 
bian matrix at that point has one zero and two purely imaginary eigenvalues. 
Consequently the Hartman-Grobman theorem does not apply. Therefore, we 
cannot draw any conclusions about the stability of the equilibrium from an 
examination of the Jacobian. 

We simplify the system by rescaling the variables by the equations 

y -> V4VW3 y, 



3 
1 

t -> . t. 

In order to take account of the equilibrium point corresponding to the point 
at "infinity" and to remove the transcendental functions, it is convenient to 
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introduce the variable u defined by 

u := e~+, 



(9) 



to obtain finally 

u = -uy, 

V = ~Hy -u(l-u), 

H=-\H 2 -\y 2 + \{l-uf. (10) 

Note that under the transformation Q, the resulted three-dimensional dy- 
namical system (fTUj) is quadratic. In view of © we have 3H 2 — \y 2 — V (0) > 
0, hence, the phase space of the system (fTUjl is the set 

E := | (u, y, H) e M 3 : H 2 - y 2 - | (1 - uf > o} . (11) 

The equilibrium points of ()10j) are: 

A: (u = 1, y = 0, H = 0) . This corresponds to the limiting state of an 
ever-expanding universe with H — > while the scalar field approaches the 
minimum of the potential and the scale factor goes to infinity. Equality in 
(|11|) which arises from the flat, k = 0, case defines a set on the boundary of S. 
We conclude that the point A which corresponds to the Minkowski solution, 
is located on this boundary. The detailed structure of this equilibrium will 
be analyzed in the next Section. 

B: [u — 0, y — ±V3/2, H = 0) . These lie outside of the phase space and, 
therefore, are unphysical. 

C: (u = 0, y = 0, H — ± a/3/2 j . In the next Section we show that only 
the point with the + sign can be approached by a trajectory starting with a 
H > 0. It corresponds to the de Sitter universe with a cosmological constant 
equal to y/V^. Regarding the stability of this equilibrium, it is easy to see 

that the Jacobian matrix of (|T0|) at q = ^0, 0, a/3/2 j has one zero and two 

negative eigenvalues. The center manifold theorem implies that there exists 
a local 2-dimensional stable manifold through q (see for example [19J). That 
means that all trajectories asymptotically approaching q as t — » oo, lie on 
a 2-dimensional invariant manifold. Since q is a nonhyperbolic fixed point, 
the topology of the flow near q is nontrivial and is characterized by a one- 
dimensional local center manifold intersecting q. In the Appendix we prove 
the following result. 
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Proposition 1 The equilibrium point q = \0, 0, \/3/2 j of MCfy is locally 
asymptotically unstable. 



4 Late time evolution 

It is easy to see that at the equilibrium point, (u — 1, y — 0, H — 0) , the 
eigenvalues of the Jacobian of (fTU|) are ±i, and therefore, we cannot infer 
about the stability of the equilibrium. Nevertheless, it is the most interesting 
case, because in all other equilibria the scalar field reaches the flat plateau, 
which is impossible if we restrict ourselves to initial values of H smaller than 



y/Voo. The study of the qualitative behavior of a dynamical system near a 

nonhyperbolic equilibrium point is difficult even in two dimensions. We find 
the so-called normal form (cf. jTH] for a brief introduction) of the system (|T0|) 
near the equilibrium point (u — 1, y — 0, H — 0) . The idea of the normal form 
theory is the following: Given a dynamical system with equilibrium point at 
the origin, x = Ax + f (x) , where A is the Jordan form of the linear part 
and f (0) = 0, perform a non-linear transformation x — » x + h (x) , where 
h (x) =0 (|x| 2 ) as |x| — > 0, such that the system becomes "as simple as 



To write the system in a form suitable for the application of the normal 
form theory, we shift the fixed point to (0, 0, 0) by setting x = u — 1 and the 
system becomes 




possible" . 



x 



-y - xy, 
x + x 2 — Hy, 

1 2 

-x y — 

2 3 



y 



H 



3 



(12) 



We now perform the non-linear transformation 



x -> x - y 2 + -Hy, 

y -> y + xy + ~Hx, 
7 

H -> H + —xy, 
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and keeping only terms up to second order we obtain the system 

1 

x = -y- -Hx, 

H = -±(x 2 + y 2 )- 1 -H 2 . (13) 

Note that the results are valid only near the origin. 

Passing to cylindrical coordinates (x — r cos 9,y = r sin 9,H = H) , we 
have 

r = — rH, 
2 ' 

= 1, 

H = -— r 2 - -H 2 . (14) 
12 3 y 1 

We note that the 9 dependence of the vector field has been eliminated, so 
that we can study the system on the r — H plane. The equation 9 = 1 means 
that the trajectory in the x — y plane spirals with angular velocity 1. It is 
convenient to rescale the variables by 

r -> 6r, H -> 3H (15) 

so that the projection of ()14j) on the r — H plane is 

3 

r = — rH, 
2 ' 

H = — r 2 — H 2 . (16) 

This system belongs to a family of systems studied in 1974 by Takens |2(Jj . 

It is easy to obtain the phase portrait of (fTB|) via numerical integration. 
However, we can analyze the qualitative behavior of the trajectories using 
theoretical arguments. Firstly, (fTB*j) is invariant under the transformation 
t — ► —t, H — > —H (which implies that all trajectories are symmetric with 
respect to the r axis) and the line r = is invariant. Secondly, the system 
(I16|) has invariant lines H = cr. To see this, write 

dH —r 2 — c 2 r 2 — 1 — c 2 _^ ±^ (17) 

dr —\ar 2 — |c 
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Taking the dot product of the vector field (— \rH, — r 2 — iJ 2 ) with the 
radial vector (r, H) ± along the line H = cr we find that it is negative for 
H > and positive for H < 0. Therefore, the direction of the flow along 
H = cr in the first quadrant is towards the origin and goes away from the 
origin in the second quadrant. Note that H is always decreasing along the 
orbits while r is decreasing in the first quadrant. Since no trajectory can cross 
the line H = cr, all trajectories starting above this line, approach the origin 
asymptotically. On any orbit starting in the first quadrant below the line 
H = cr, H becomes zero at some time and the trajectory crosses vertically 
the r— axis. Once the trajectory enters the second quadrant, r increases and 
H decreases. The phase portrait is shown in Figure 1. At first sight, it seems 




Figure 1: The phase portrait of (fTBj) . The invariant trajectories H = ±^/2r 
are shown with thicker lines. For expanding models only trajectories above 
the line H = y/Qr belong to the phase space of the system. 

probable that an initially expanding universe may recollapse. However, the 
phase space of the dynamical system (|TB|) is not the whole r — H plane, 
because of the constraint (jllj) , which in terms of the variables (II 5|) becomes 

H 2 > 6r 2 . (18) 
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Therefore, for an expanding universe we should consider only trajectories 
starting above the line H = a/oV and according to the previous discussion 
all these trajectories asymptotically approach the origin. 

We now turn to the relation of the dynamics in the r — H plane to 
the full three-dimensional system (jl4j) . or the equivalent (|T3|) in Cartesian 
coordinates. Any trajectory spirals clockwise in the x — y plane while both H 
and x 2 + y 2 are decreasing. In physical terms this means that the scalar field 
oscillates around the minimum of the potential with a decreasing amplitude, 
H is always decreasing and, in view of ©, the curvature decreases. A typical 
trajectory of (fT3j) is shown in Figure 2. 




Figure 2: Trajectories of the full three-dimensional system spiral approaching 
the origin. 



5 Comment on Flat FRW models 

Consider a flat FRW model containing a barotropic fluid with an equation of 
state p = (7 — where < 7 < 2, and a scalar field having the potential 
((7j). Then the system ©-(01) reduces again to a four-dimensional dynamical 
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system, namely 



subject to the constraint 



y 
p 

H 



y, 

-3Hy - V (0) , 
-3 1P H, 

1 2 7 
~2 V ~ 2 P ' 



(19) 



3H Z 



v(<j>). 



(20) 



In contrast to (0i), the forth of ()19|) implies that H is always decreasing. If 
we use the constraint (J2*U|) to eliminate p, the evolution equations become 



= y, 

y = -3Hy-V'{(f>) 
H = 



-—H 2 
2 



4 y 2 



(21) 



We see that the structure of (}2~Tj) describing a flat FRW model with a perfect 
fluid and a scalar field, apart from being parameter dependent, has a striking 
similarity to the dynamical system (jHJ) for the vacuum positively curved FRW 
with a scalar field. Proceeding as in Section IV, we end up with the following 
system in cylindrical coordinates 



r 

H 



--rH, 

2 ' 



5 7~ 4 r 2 _ X H 2 



(22) 



Although the system (|22|) depends only on one parameter, it is convenient 
to rescale the variables by 



with 



Ar, H — > fiH 



16 



(23) 



7 



7 (57 - 4) 
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so that the projection of (|22|) on the r — H plane is 



r = rH, 

7 

H = br 2 -H 2 , (24) 



where 



b = +1 for 7 > 4/5, 6 = -1 for 7 < 4/5. 
Note that ()24)1 has a first integral, viz. 

, (r ,„ ) = _i_ r -^(_^__^). ( 25 ) 

In fact, it is straightforward to verify that (dl/dr) r + (dl/dH) H = along 
the solution curves of ()24|) . The level curves of / are the trajectories of the 
system (cf. Figure 3). 

Invariant lines H = cr exist for certain values of the parameter 7. We 
find (cf. (HZD) 



±4 ' 1-1/7- 

Case I. For b = +1, invariant lines exists if 7 > 1. We find that the 
direction of the flow along H = cr in the first quadrant is towards the origin 
and goes away from the origin in the second quadrant. Note that in the first 
quadrant r is decreasing along the orbits and that H vanishes along the line 
H = r, which lies below the invariant line H = cr. It can be shown that in 
the first quadrant a level curve of / (r, H) may intersect the line H = r only 
once (it is sufficient to consider the level curve passing through an arbitrary 
point (ri,0) and compute the r coordinate at the intersection with the line 
H = r). We conclude that once a trajectory crosses the line H = r, it is 
trapped between the lines H = r and H = cr and, since f < 0, it approaches 
the origin asymptotically. 

Case II. b = +1 and 4/5 < 7 < 1. There are no invariant lines. Similar 
arguments as in case I, yield the phase portrait shown in Fig. 3. 

Case III. 7 < 4/5 (=^> b = —1) . The analysis is exactly the same as in 
Section IV. In all cases we must remember that the phase space of the dy- 
namical system (J24j) is not the whole r — H plane, because of the constraint 
(II which in terms of the variables (|23|) reads 

H2 > ^tV- < 26 > 
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Figure 3: The phase portrait of (|2"4*j) for 1 < 7 < 2, 4/5 < 7 < 1 and 
< 7 < 4/5 respect ivelly. 

Therefore, for an expanding universe, we should consider only trajectories 
starting above the line H = A/67/ [5~7 — 4|r, which in case III lies always 
above the line H = cr and according to the previous discussion all these 
trajectories asymptotically approach the origin. 

We conclude that, in the conformal frame of the R + aR 2 theory, an ini- 
tially expanding flat universe with a barotropic fluid as matter source remains 
ever-expanding and eventually the quadratic curvature corrections become 
negligible. This result, established by stability analysis, is in accordance 
with the general properties of all flat FRW models with a scalar field having 
a potential with a unique zero minimum [T7j . 

6 Discussion 

We have analyzed the qualitative behavior of a positively curved FRW model 
containing a scalar field with the potential (J7J). This model is conformally 
equivalent to the positively curved FRW spacetime in the simplest higher 
order gravity theory, namely the R + aR 2 theory. We have shown that 
an initially expanding closed universe avoids recollapse provided that the 
initial value of H is less than y/Voo- This result should be compared to a 
counterexample of the closed universe recollapse conjecture (cf. [21] where 
it is shown that initially expanding vacuum diagonal Bianchi IX models in 
purely quadratic gravity are ever-expanding). It should be of interest to 
investigate if a closed FRW universe filled with ordinary matter satisfying 
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the usual energy conditions and a scalar field with the potential (J7J) can avoid 
recollapse. This is equivalent to the analysis of the qualitative behavior of 
the full five-dimensional system (j2J)-©- A partial answer to this question for 
an arbitrary non-negative potential having a unique minimum V (0) = is 
given in [T7] . 
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mous referee for useful suggestions, especially for pointing me out the remarks 
regarding the equilibrium A after eq. 



A Proof of Proposition 1 

In order to determine the local center manifold of (jlO)) at q, we have to trans- 
form the system into a form suitable for the application of the center manifold 
theorem. The procedure is fairly systematic and will be accomplished in the 
following steps. 

1. The Jacobian of (fTUj) at q = ^0,0, a/3/2^ has eigenvalues 0, —a/2/3 

and —a/3/2 with corresponding eigenvectors ^— a/2/3, 2/3, 1 j , (0,0,1) 

and (0, 1, 0) . Let T be the matrix having as columns these eigenvectors. 
We shift the fixed point to (0, 0, 0) by setting H = H- y/3/2 and write ((TOD 
in vector notation as 

z = Az + F(z), (A.l) 

where A is the linear part of the vector field and F (0) = 0. 

2. Using the matrix T which transforms the linear part of the vector field 
into Jordan canonical form, we define new variables, (x, 2/1,2/2) = x, by the 
equations 

2 

y = ^ x + 

H = x + y 1 , 

or in vector notation z = Tx, so that (jA.l|) becomes 

x = T^ATx + T _1 F (Tx) . 
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Denoting the canonical form of A by B we finally obtain the system 



where f fx) 



x = Bx + f(x), 
T _1 F (Tx) . In components system (|A.2|) is 



(A.2) 





/o 


X \ 





*■ = 


2/2 / 


1° 












+ 







2 9 

-far - xy 2 

10^2 _ 12 _ 

27 X 3»1 
4 9 2 

-x z - ~ xyi - 



\ 



2 / 



3 2/2 
3^2 




3. The system 



9"^ 3 

3j) is written in diagonal form 

x = Cx + f (x, y) 
y = ^Py + g (x, y) , 



- 2/12/2 



(A.3) 



(A.4) 



where (x, y) G R x IR 2 , C is the zero lxl matrix, P is a square matrix 
with negative eigenvalues and /, g vanish at and have vanishing derivatives 
at 0. The center manifold theorem asserts that there exists a 1-dimensional 
invariant local center manifold W c (0) of (jA.4|) tangent to the center subspace 
(the y = space) at 0. Moreover, W c (0) can be represented as 



W c (0) = {(x,y)G 



x 



: y = h(x), \x\ < 5} ; h (0) = 0, Dh (0) = 0, 



for 5 sufficiently small (cf. ^J], P- 155). The restriction of (jA.4j) to the 
center manifold is 

x = Cx + f(x,h(x)). (A.5) 

According to Theorem 3.2.2 in [22j, if the origin x = of (jA.5|) is stable (resp. 
unstable) then the origin of (jA.4|) is also stable (resp. unstable). Therefore, 
we have to find the local center manifold, i.e., the problem reduces to the 
computation of h (x) . 

4. Substituting y = h (x) in the second component of (jA.4|) and using the 
chain rule, y = Dh (x) x, one can show that the function h (x) that defines 
the local center manifold satisfies 



Dh (x) [Cx + / (x, h (x))} - Ph (x) - g (x, h (x)) = 0. 



(A.6) 
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This condition allows for an approximation of h (x) by a Taylor series at 
x = 0. Since h (0) = 0, Dh (0) = 0, it is obvious that h (x) commences with 
quadratic terms. We substitute 

i / \ _ ( hi (x) \ _ / aix 2 + a 2 x 3 + O (x 4 ) \ 
{X) ~ : \h 2 {x) J ~ V hx 2 + b 2 x 3 + O (x 4 ) J 

into (jA.6|) and set the coefficients of like powers of x equal to zero to find the 
unknowns a 1; 6 1; .... 

5. Since yi is absent from the first of (jA.3|) . we give only the result for 

h 2 (x) . We find h = f^f, b 2 = ±. Therefore, (lA31) yields 

x = -^ 2 -^x 3 -^x 4 + 0(x 5 ). (A.7) 

It is obvious that the origin x = of (|A.7|) is asymptotically unstable (saddle 
point). The theorem mentioned after (|A.5|) implies that the origin x = of 
the full three-dimensional system is unstable. This completes the proof. 
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